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Nanotubular materials inspired by crystalline diborides such as AIB2 are proposed. The atomic 
structure, in particular the basic chemical question of where to put Al atoms in order to stabi- 
lize nanotubular Al-B systems, is investigated using density-functional calculations for prototype 
systems. The optimized tubular prototypes are found to be competitive in energy with their bulk 
crystalline counterparts. All of the tubular Al-B systems investigated are calculated to be metallic. 

PACS numbers: 61.46.+W, 73.22.-f, 73.63.Fg, 81.07.Dc 
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I. INTRODUCTION 

Soon after the discovery of carbon nanotubesJH it was 
recognized that these compounds can have electronic 
properties ranging from semiconducting to metallic, de- 
pending on the details of their atomic structure (see chap- 
ter 19 in Ref. |^). This has led to much speculation about 
nanotubular systems becoming the materials of choice 
for future miniaturizations of electronic devices towards 
the nanodomain.EI Some major challenges are the con- 
trol of the growth mechanism in order to obtain carbon 
nanotubes with specific electronic properties, and the de- 
velopment of techniques to wire together semiconducting 
and metallic nanotubes in a controlled manner. One ap- 
proach for connecting carbon nanotube devices could be 
the use of other metallic nanotubes that function as wires 
on the same length scale. 

In an earlier publication^ we discussed the possible 
existence of metallic boron nanotubes that could serve 
such a purpose. Boron nanotubes can be constructed 
from the so-called Aufbau principle for the formation of 
stable boron compounds formulated in Ref. ^. This 
Aufbau principle identifies a small set of simple struc- 
tural features characteristic of the most stable isomers of 
pure boron clusters and proposes that highly stable boron 
clusters, surfaces, and networks can be constructed using 
only these structural elements. Crystalline a-boron has 
long been the paradigm for boron bonding. It is com- 
prised of Bia_ icosahedra with six-fold inverse-umbrella 
coordination)!] a configuration that obeys the Aufbau 
principle. The coordination in the proposed boron nan- 
otubes is very different - hexagonal pyramidal - but also 
follows the Aufbau principle. In a recent study of clus- 
ters of 96 boron atoms J3 we showed that the formation 
of large quasi-planar and tubular systems with six-fold 
pyramidal bonding is favored over the formation of clus- 
ter aggregates typical for crystalline a-boron. We also 
found that the strain energy of boron nanotubes rela- 



tive to the corresponding sheets is intermediate, between 
those of the carbon and boron-nitride systems Jj another 
favorable condition for the formation of boron nanotubes. 



In this paper, the large family of layered metal di- 
borides is used as a starting point for investigating metal- 
boron nanotubes. Diborides of the AIB2 lattice type have 
graphene-like layers of boron atoms, with metal atoms 
located halfway between the boron layers. These struc- 
tures are similar to the pure boron surfaces discussed in 
Ref. ||. In fact, in Ref. |, the AIB2 system was used to 
help elucidate the quasi-planar configurations found for 
pure boron. Here we extend the analogy to the tubu- 
lar systems. Density functional calculations are used to 
explore the geometry, stability, and electronic structure 
of some prototype Al-B nanotube materials. The pro- 
totype nanotubular systems presented here are found to 
be energetically competitive with the corresponding bulk 
counterparts, suggesting that synthesis of such nanotubes 
might be possible. 



The remainder of this paper is organized as follows. 
Sec. II is devoted to the structural properties of the di- 
boride systems. We analyze the A1B 2 lattice type, de- 
scribe some of the basic physical properties of diborides, 
and illustrate the geometric construction of metal-boron 
nanotubes starting from metal-diboride systems. In 
Sec. Ill, the technical details of the ab initio calcula- 
tions are given. In Sec. IV, we present some prototypical 
Al-B nanotubular structures that are calculated to be 
energetically favorable, and we discuss the bonding and 
electronic structure of these systems. Conclusions and 
remarks about future work are given in Sec. V. The 
present study proposes a new class of nanotubular mate- 
rials and opens the door to a range of interesting prob- 
lems for future studies. 



2 



II. CONSTRUCTION OF NANOTUBES FROM 
CRYSTALLINE ALB 2 

The typical diboride structure is the AIB2 lattice, 
which is composed of two parallel systems of flat layers.cl 
One layer contains boron atoms arranged in a honey- 
comb lattice like a graphene sheet (see Fig.|l|(a)). The Al 
atoms form a triangular lattice that is located halfway 
between the graphene-like sheets, positioned so that the 
Al atoms project down onto the centers of the boron 
hexagons (see Fig.^(b)). The space group is P6/mmm, 
and the primitive cell is hexagonal with one formula unit 
per cell. Topically, the in-plane lattice constant is around 
a = 3.1 AB This corresponds to an in-plane B-B distance 
of about 1.8 A, which is similar to the B-B bondlcngth 
in a-boron. The interplanar B-B distance ia. significantly 
larger, with typical c/a ratios of about 1.1 B 

The family of diborides of the AIB2 type comprises 
materials with a wide variety of metal atoms substitut- 
ing for the_Al position. There are at least 22 members of 
the family,Q ranging from MnB2 with an effective metal- 
lic radius of 1.29 A to GdB2 with an effective metallic 
radius of 1.81 A. Simple sp elements, as well as d- and /- 
electron atoms can be accommodated. The ability of the 
AIB2 lattice to incorporate metals that differ so much in 
size and electronic structure is remarkable, and it raises 
the possibility of forming alloys between members of this 
family. 

The construction of nanotubular structures from the 
AIB2 lattice type is shown in Fig. |l|. Ignoring the ampli- 
tudes of the Al atoms for the moment (Figjl](b)), we see 
that A1B 2 nanotubes can be generated by the same cut- 
and-paste procedure used to generate carbon nanotubes 
from graphene sheets.0 Nanotubes may be classified by a 
pair of integers ( N, M) , that specify a cut along a vector 
Nhi + Mb2, where bi and b2 are primitive lattice vec- 
tors for the honeycomb lattice. This first cut is followed 
by a second cut perpendicular to this direction that con- 
tinues until another atom of the honeycomb lattice is hit. 
Figure |l] shows the resulting sheet for a (6,6) type of nan- 
otube in both systems. Rolling up a general sheet and 
pasting it along the perpendicular direction from both 
ends of the (N,M) cut vector generates the basic unit 
cell of an (N, M) nanotube. For the purpose of ab initio 
simulations, it is common practice to place these tubes 
side by side on a hexagonal superlatticepca thereby ap- 
proximating the bundles in which these compounds form 
in Nature. 

It is interesting to note that if B is substituted for Al 
in Fig. ^(b), the structure becomes a puckered layer of 
six-fold coordinated boron. If the construction described 
above is applied, we obtain the pure boron nanotubes 
described in Ref. ^, where we also introduced a different 
but equivalent construction procedure labelled by a pair 
of integers that describes cuts along the unit vectors of a 
triangular lattice. We conjectured that the best way to 
put boron nanotubes on a hexagonal superlattice would 
be to arrange them such that pairs of boron atoms meet 



along the directions of the underlying hexagonal lattice.cl 
Such an arrangement leads to a type of chemical bonding 
between the tubes that is analogous to what is found 
for planar forms of pure boron. In the planar systems, 
pairwisc bonds between 7-fold coordinated boron sticking 
out of puckered 6-fold coordinated layers of boron were 
found to be the most favorable arrangement for binding 
single sheets together to form layered systems. Ej 

In earlier studiesjj we found that tubular boron sys- 
tems could stand alone, i.e., like carbon nanotubes, they 
like to bundle in the form of a hexagonal superlattice, but 
they do not have to. For the diboride system, the situa- 
tion may be different. Studies of small boron clusters in- 
dicate that neither the honeycomb lattice nor the related 
carbon-type nanotube structure is stable for boron.El In 
fact, the local atomic structure of the diboride systems 
suggests that metal atoms are necessary to stabilize hon- 
eycomb boron layers or the tubes derived from such lay- 
ers. The metal atoms themselves have to bind to a suffi- 
cient number of boron atoms (12 in the case of the AIB2 
lattice). It is not clear a priori whether one could sim- 
ply bind metal atoms inside or outside of a carbon-type 
boron tube to stabilize a single tube, or whether those 
tubes are more likely to come in larger bundles of carbon- 
type boron nanotubes with metal atoms both inside and 
between the tubes, serving in part to bind the tubes to- 
gether. For Al-B tubes, our ab initio calculations show 
that the latter arrangement is preferred, as discussed in 
Sec.IV. 



III. TECHNICAL DETAILS 

Calculations were carried out,.-jJsing the VASP 
package,tM a density functionalE3tL3 based ab initio 
total-energy code usingnplane-wave basis sets and ultra- 
soft pseudopotentialsJlat-a When performing structural 
optimizations, we let all degrees of freedom relax, i.e., 
the complete set of atomic positions as well as the pa- 
rameters of the unit cell. The optimal configurations as 
well as the electronic structure for each relaxation step 
were determined using preconditioned conjugate gradient 
procedures .EjEj 

The kinetic-energy cutoffs used for the plane-wave ex- 
pansion of electronic wavefunctions were 358.2 eV for C 
(which we use as a reference system), 161.5 eV for Al, and 
321.4 eV for B and all of the Al-B systems. The Brillouin 
zone was sampled on grids of (13x13x13) k-points in the 
cases of crystalline AIB2 and fee Al, (8x8x8) k-points in 
the case of graphite, and (5x5x5) k-points in the case of 
all tubular systems with the exception of the (6,6) carbon 
nanotube for which (4x4x4) k-points were used. These 
meshes yield meV accuracy in the binding energies, and 
their differing sizes reflect differences in unit-cell sizes. 
Care was taken to include a sufficient number of bands 
above the Fermi level to avoid numerical problems in the 
iterative diagonalization scheme. 

Although the calculations were well converged with re- 
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spect to the number of plane waves and k-points, the re- 
sults presented in Sec. IV have to be interpreted with 
care. In particular, the ultrasoft pseudopotential for 
boron seems to overbind a-boron by about 0.7 eV/atom 
as compared to results for the same system obtained 
with harder norm-conserving pseudopotentials.E3t3 On 
the other hand, the cohesive energy we calculate for the 
diboride AIB2 is in excellent agreement with receat all- 
electron LDA calculations on diboride systems .Ej This 
suggests that the ultrasoft pseudopotential for boron has 
difficulty describing the inverse-umbrella bonding in a- 
boron, but it nevertheless gives an accurate description 
for nearly planar configurations of boron like the ones 
found in the diborides, quasi-planar boron sheets, and 
their related nanotubes. 

Another potential problem is the fact that the conju- 
gate gradient optimization does not guarantee the detec- 
tion of a global minimum. Given the reported pathol- 
ogy of the energy hypersurfaces of boron clusters,u it is 
likely that for the large systems discussed below, we relax 
into local minima rather than global minima. However, 
as discussed in Ref. [|, it is also likely that the relaxed 
structures capture the right chemistry, and all we have 
to worry about is the fact that the cohesive energies cor- 
responding to these local minima will be slightly higher 
than those of the global minimum. 



IV. RESULTS AND DISCUSSION 

Unlike pure boron and pure carbon nanotubes, the Al- 
B nanotubes we investigated do not appear to be stable as 
stand-alone tubes. We were unable to converge structural 
optimizations of isolated tubes constructed by rolling up 
puckered sheets of AIB2. Furthermore, in structural op- 
timizations of Al-B nanotubes on periodic superlattices, 
the tubes always ended up packed closely together, even 
if they started out well separated. Hence, for the re- 
mainder of this paper, we focus on results for Al-B tubes 
arranged on a hexagonal superlattice. For comparison, 
parallel calculations were performed for pure boron and 
pure carbon nanotubes. In each case, the (6,6) nanotube 
was used as a prototype. 

The first step was to generate reference data for atomic 
structure and cohesive energies. This data is listed in 
the lower part of Table |. Fee Al and a-boron serve as 
basic reference structures for the Al-B systems. Phase- 
separated AIB2 (into fee Al and a-boron) is calculated 
to have a cohesive energy of 6.38 eV/atom. As can be 
seen in Table |, crystalline AIB2 is calculated to be only 
slightly stable against phase separation. However, the va- 
lidity of this result is questionable in light of the known 
stability of crystalline AIB2.Q As we did a full relaxation 
on all systems, with high plane-wave cutoffs and many 
k-points, we attribute the problem to the overbinding of 
a-boron discussed in Sec. III. As already mentioned, 
the cohesive energy— pf AIB2 obtained in a recent all- 
electron calculationEj is in excellent agreement with our 



value listed in Table |. 

The same considerations apply when we examine the 
stability of the (6,6) boron nanotube. The results in Ta- 
ble [l| indicate that the (6,6) boron nanotube is unstable 
by 0.5 eV/atom compared to crystalline a-boron. This 
is in contrast to our recent all-electron ab initio calcula- 
tions, which found that for Bgg clusters, the nanotubular 
isomers challenge or exceed the stability of cluster ag- 
gregates found in a-boron.0 Again, we attribute this dis- 
crepancy to the overbinding of a-boron by the ultrasoft 
pseudopotential used. 

The situation looks more promising when we analyze 
the chemistry of the boron nanotubes. The relaxed struc- 
ture for the (6,6) boron nanotube is shown in Fig. |2|(b). 
Note that the representation of the tubes in Fig. g is per- 
haps unusual in that no tube is shown in its entirety. In- 
stead, this representation is used to emphasize the chem- 
ical bonding between nanotubes sitting on the hexagonal 
superlattice. Prominent boron to boron bonding along 
the direction of the lattice is evident in Fig. |||(b) . Fur- 
ther, there is a slight faceting due to undulations from 
the basic hexagonal lattice of the quasiplanar reference 
structure. Both of these features are in perfect agree- 
ment with the structure and the chemical bonding found 
in boron sheets.B Therefore, it is reasonable to take the 
relaxed (6,6) boron nanotube as a good starting structure 
when looking for nanotubes of the Al-B type. 

Table | includes an additional set of benchmarks. 
By comparing graphite and the (6,6) carbon nanotube 
(Fig. ||(a)), we get a realistic measure for the energetic 
difference between comparable carbon nanotubes and 
their layered reference structure. We calculate the (6,6) 
carbon nanotube to be about 0.14 eV/atom higher in en- 
ergy than graphite. Tubes with larger radii would have 
even smaller instabilities .□ 

Among the many variants of Al-B nanotubes that we 
explored, we focus on the two shown in Figs, ^(c) and (d). 
Our approach was to start with the relaxed B nanotube 
and substitute Al for B. The structure shown in Fig. @(c) 
has composition A1 3 B 30 and was generated from the (6,6) 
nanotube by removing pairs of boron atoms located along 
the directions of the hexagonal lattice, and substituting 
them by single Al atoms. This can happen at three dif- 
ferent places within the unit cell, and the Al atom then 
sees a local atomic environment similar to AIB2. Indeed, 
after full relaxation, the system tends to recreate these 
local environments, and the tube as a whole becomes less 
puckered (compare Fig. ||(c) to the puckered B tube in 
Fig. ^(b) and the unpuckered C tube in Fig. ||(a)). How- 
ever, the energetic gain is not substantial: the cohesive 
energy of phase-separated AI3B30 is about 7.2 eV/atom, 
so there is only a small improvement in stability as com- 
pared to the (6,6) boron tube. 

Figure ||(d) shows an Al-B nanotube system with sto- 
ichiomctry AIB2. This configuration is calculated to be 
unstable by only 0.23 eV/atom compared to crystalline 
AIB2. While this energy difference is not quite as small 
as that between the (6,6) carbon tube and graphite, it is 
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close. The structure was generated from the (6,6) boron 
nanotube by substituting Al for boron at the six positions 
located along the directions of the hexagonal superlat- 
tice, which are on the outside layer of the puckered tube, 
and at six positions on the inside layer of the tube. In 
the optimized structure, the tubes are rotated from the 
orientations they adopt in the superlattice of pure boron 
nanotubes in such a way that Al atoms sticking out from 
adjacent tubes do not face each other, but rather sit be- 
tween boron layers. Such an arrangement locks the tubes 
together in a gear-like fashion, and creates a local envi- 
ronment similar to that in crystalline AIB2. The inclu- 
sion of Al atoms inside the tube is found to be important 
for stabilizing this structure. 

The electronic densities of states (DOS) calculated for 
crystalline AIB2 as well as the (6,6) B, AI3B30, and AIB2 
tubes are plotted in Fig. 0. All of the tubes are metallic. 
This result was predicted in Ref. || based on qualitathze. 
arguments related to hexagonal tight-binding models,E£l 
as well as on static ab initio calculations. The present 
results show that these predictions even hold after a com- 
plete structural relaxation. The overall shape of the DOS 
for crystalline and tubular AIB2 are similar due to sim- 
ilarities in local atomic environments. The DOS of the 
AIB2 nanotube, however, has much more fine structure, 
which can be understood in terms of a backfolding of the 
bands of crystalline AIB2 into the reduced Brillouin zone 
of the tubular lattice. The DOS of AI3B30 plotted in 
Fig. |^(c) on the other hand shows some similarity to the 
DOS of the (6,6) boron nanotube depicted in Fig. ||(b), 
with some distortion due to the modest content of Al 
atoms. 



V. CONCLUSIONS 

We have explored the plausibility of metal-boron nan- 
otubular systems derived from diboride materials. Our 



data suggests that it may be possible to synthesize such 
systems. We predict that metal-boron nanotubular ma- 
terials will likely form in bundles built from rolled up 
graphene-like sheets of boron, with metal atoms sitting 
inside and outside the tubes, preserving as much as pos- 
sible the local atomic environments found in the AIB2 
lattice. Ab initio calculations of their basic electronic 
properties show that such tubular metal-boron systems 
are metallic. 

There are numerous issues of interest for future studies. 
First, a comparison of the properties expected for nan- 
otubes derived from different diboride compounds - in- 
cluding an assessment of which, if any, form stand-alone 
tubes - would be of interest. Continuing along those 
lines, it may be possible to tune some of the properties of 
the metal-boron tubes by forming alloys. Alternatively, 
substitution of boron atoms by carbon is likely to have a 
serious impact on the electronic structure of the resulting; 
tubular materials, bringing them closer to BN materials.El 
Also, the recent discovery of superconductivity in MgB2 
raises the question of whether some of the._diboride- 
derived tubular materials might superconductEi Finally, 
and most importantly, we look forward to cooperating 
with experimental groups that may be able to synthesize 
tubular metal-boron materials. 
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TABLE I: Number of atoms per unit cell n, and cohesive en- 
ergies E co h, for pure and mixed boron nanotubes and several 
reference structures. 



System 


n 


E coh [eV/atom] 


(6,6) boron nanotube 


36 


7.001 


AI3B30 nanotube 


33 


6.884 


AIB2 nanotube 


36 


6.227 


a-boron 


12 


7.509 


fee Al 


1 


4.119 


A1B 2 


3 


6.451 


graphite 


4 


10.156 


(6,6) carbon nanotube 


24 


10.019 
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(a) Generating a (6,6) carbon nanotube 














b2* ° * 







(b) Generating a template (6,6) A1B 2 nanotube 



FIG. 1: Cut-and-paste construction for (a) (6,6) carbon nan- 
otubes and (b) template (6,6) AIB2 nanotubes. In (a) circles 
[o] denote C sites; in (b) circles [o] denote B sites and crosses 
[+] denote the projection of Al sites onto the B plane. In each 
(6,6) tube is generated by rolling up the shaded strip, 
the edges of which are defined by the vector 6bi + 6b 2 and 
the shortest lattice vector perpendicular to this direction. 



(a) (6,6) carbon nanotubo 



(b) (6,6) boron nanotubo 




(c) AI3B30 nanotube (d) AIB2 nanotubo 



FIG. 2: Comparison of various optimized tubular structures. Each panel shows one unit cell of the hexagonal superlattice on 
which tubes are arranged. In (c) and (d), the light atoms are B and the dark atoms are Al. Units of the surrounding boxes are 
given in A. Note that even after relaxation, the atoms remain stacked in two layers only. 
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FIG. 3: Densities of state (DOS) in states/eV/cell vs. energy in eV for crystalline AIB2 and various optimized tubular systems. 
All systems are metallic. Similarities in the DOS arise from similarities in local atomic environments. 



